This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 



BEST AVAILABLE IMAGES 



Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not limited to): 



BLACK BORDERS 

TEXT CUT OFF AT TOP, BOTTOM OR SIDES 
FADED TEXT 
ILLEGIBLE TEXT 
SKEWED/SLANTED IMAGES 
COLORED PHOTOS 

BLACK OR VERY BLACK AND WHITE DARK PHOTOS 
GRAY SCALE DOCUMElNfTS 



IMAGES ARE BEST AVAILABLE COPY. 

( 
) 

As rescanning documents will not correct images, 
please do not report the images to the 
Image Problem Mailbox. 



THIS PAGE BUNK (uspto) 



612 .40181X00 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE ^ 



Applicant (s) : DURET, et al . 



<30 0\ 



Serial No.: (Not yet assigned) So go i 



Filed: May 10, 2001 

Title: AUTOMATIC PIPE GRIDDING METHOD ALLOWING 

IMPLEMENTATION OR FLOW MODELLING CODES 

Group: (Not yet assigned) 



LETTER CLAIMING RIGHT OF PRIORITY 

Honorable Commissioner of „ June 19, 2001 

Patents and Trademarks 
Washington, D.C. 20231 

Sir : 

Under the provisions of 3 5 USC 119 and 3 7 CFR 1.55, the 
applicant (s) hereby claim (s) the right of priority based on 
French Patent Application No.(s) 00/08200, filed June 23, 
2000 . 

A certified copy of said French Application is attached. 

Respectfully submitted, 

ANTONELLf\ TERRY, STOUT & KRAUS , LLP 



Donald E. Stout 
Registration No. 26,422 



DES/alb 
Attachment 
(703) 312-6600 



THIS PAGE BLANK (uspto) 



I INSTITUT 
NATIONAL DE 
LA PROPRIETE 
INDUSTRIELLE 



REVET D ' I N VE N T I O 



• CERTIFIGAT D'UTILITE - CERTIFICAT D'ADDITION 



COPIE OFFIGIELLE 



Le Directeur general de NnstiM la propriete 

industrielle certifie que le document ci-annexe est la copie 
certifiee conforme d'une dennande de titre de propriete 
industrielle deposee a I'lnstitut. 



P , p ■ , 2 8 MAI 2001 

Fait a Pans, le - 



Pour le Directeur general de I'lnstitut 
national de la propriete industrielle 
Le Chef du Departement des brevets 




Martine PLANCHE 



INSTITUT 
NATIONAL DE 
LA PROPRIETE 
INDUSTRIELLE 



SIEGE 

26 bis, rue de Saint Petersbour^ 
75800 PARIS cedex 08 
Telephone : 01 53 04 53 04 
Telecopte : 01 42 93 59 30 
http;//www.tnpi-fr 



ETABLISSEMENT PUBLIC NATIONAL CREE PAR LA LOt N" 51-444 DU 19 AVRIL 1951 



THIS PAGE BLANK (uspto) 




E 



3 mariTUT 

NAflONAL DB 
LAPDOPOIBTO 
INDUOiailLLI 



26 bis, rue de Saint P^tersbourg 
75800 Paris Cedex 08 

Telephone : 01 53 04 53 04 T6l6copie : 01 42 94 86 54 



REMISE DES PItCES 
DATE 

LIEU 



CEOTOFDCM D'yTDLDTE 11354-01 

Code de la propri6t6 intellectuelle - Livre VI 

REQUITE EN DtUVRANCE 1/2 

Cet imprim6 est a remplir lisiblement a I'encre noire db sao w /26os99 



D'ENREGISTREMENT 
NATIONAL ATTRIBUe PAR LINPI 

DATEDED£P6TArn^Uj Jjjjj^ 2000 
PAR L'lNPI 



0008^00 



Vos references pour ce dossier 

(factdtatif) JC/CLN 



[a] NOM ET ADRESSE DU DEMANDEUR OU DU MANDATAIRE 
A QUI LA CORRESPONDANCE DOIT EIRE ADRESSEE 

"iNSTITUT FRANCAJS DU PETROLE 
1 et 4 avenue de Bois Prdau 
92852 Rueil Malmaison cedex 



Confirmation d'un d6p6t par tdi^copie Fl N** attribu6 par I'lNPI a la telecopie 



[U NATURE DE LA DEiViANDE 



Demande de brevet 



Demande de certificat d'utilite 



Cochez i'une des 4 cases suivantes 



0 



□ 



Demande divisionnaire 

Demande de brevet initiale 
OU demande de certificat d'utilite initiale 



□ 



Date 
Date 



1 L 



I L 



Transformation d'une demande de 

b revet eu ropeen Demande de brevet initiale 



□ 



Date L 



_L 



] TITRE DE L'lWVEWTlOW (200 caracteres ou espaces maximum) 

METHODE DE MAE^LAGE AUTOMATIQUE DE CONDUITES PERMETTANT L' IMPLEMENTATION DE CODES DE 
MODELISATION DE FLUX 



DECLARATION DE PRIORITE 
OU REQUETE DU BI^N^FICE DE 
LA DATE DE DEPOT D'UNE 
DEMANDE ANT^RIEURE FRANQAISE 



Pays OU organisation 
Date L_Z l—^ 



Pays OU organisation 
Date I L 



Pays OU organisation 

Date I / / I 



I I S'il y a d'autres priorit^s, cochez la case et utilisez rimprime «Suiteu 



DEMANDEUR 



I I SMI y a d'autres demandeurs, cochez la case et utilisez rimprime (fSuite» 



Nom OU denomination sociale 



INSTITUT FRANCAIS DU PETROLE 



Prenoms 



Forme juridique 



Organisme Professionnel 



N° SIREN 



Code APE-NAF 



Adresse 



Rue 



1 & 4 avenue de Bois Preau 



Code postal et ville 



92852 



Rueil Malmaison cedex 



Pays 



France 



Nationalite 



Frangaise 



N° de telephone (facultatif) 



01 47 52 60 00 



N° de telecopie (facultatif) 



0147 50 70 03 



Adresse electronique (facultatif) 



BREVET AhVENTION fr. 
iN^I CERTIFICAT D'UTILITE 



■ INSTITVT 
HATIONAL OB 



REQU£TE EN DEUVRANCE 2/2 



REMISE OES PIECES 
DATE 

LIEU 



R6serv6 ^ riNPt 



N** D'ENREGISTREMENT 
NATIONAL ATTRIBUt PAR L'lNPI 



OOOS &0 



pa 5<Q W /260e99 



Vos references pour ce dossier : 

(facultatij) 


JC/CLN 


□ MANDATAIRE 




Nom 


ELMALEH 


Prenom 




Cabinet ou Societe 


DVSTTTUT FRANCAIS DU PETROLE 


N °de pouvoir permanent et/ou 
de lien contractuel 




Adresse 


Rue 


1 & 4 avenue de Bois Preau 


Code postal et ville 


92852 Rueil Malmaison cedex 


N° de telephone (facultatij) 


01 47 52 60 00 


N° de telecopie (facultatij) 


01 47 52 70 03 


Adresse electron i que (facultatij) 




B INVENTEUR (S) 




Les inventeurs sont les demandeurs 


nOui 

1 X 1 Non Dans ce cas fournir una designation d'inventeur(s) s6par6e 


Q RAPPORT DE RECHERCHE 


Uniquement pour une demande de brevet (y compris division et transformation) 


Etablissement immediat 
ou etablissement differe 


s 

□ ^. 


Paiement echelonne de la redevance 


Paiement en deux versements, uniquement pour les personnes physiques 

□Oui 
ITlNon 


El REDUCTION DU TAUX 
DES REDEVANCES 


Uniquement pour les personnes physiques 

CDRequise pour ta premiere fois pour cette invention (joindre un avis de non-imposition) 

C^Requise anterieurement a ce dep5t (joindre une copie de la decision d'admission 
pour cette invention ou indiquersa reference) : 




Si vous avez utilise rimprimd uSuite»r 
indiquez le nombre de pages jointes 






EE SIGNATURE DU DEMANDEUR 

OU DU MANDATAIRE ^^'^"O 
(Nom et qualite du signataire) y^^^'^^^^^ 
Alfred ELMALEH, C y>^^^ 
Chef du Departement Brevets ^^^^ 
N° 422-5/PP.253 


VISA DE LA PREFECTURE 
OU DE L'lNPI 



La loi n°78-17 du 6 Janvier 1978 relative a I'informatique, aux fichiers et aux libertes s'applique aux reponses faites a ce formulalre. 
Elle garantit un droit d'acces et de rectification pour les donn6es vous concernant aupres de I'lNPI. 



TSpi 



■ INailTUT 
MATIONAL Dl 
LA PROPRirrC 
IHOUSTB'JStLI 



BREVElVlNVENTION 

CERTIFICAT D'UTILITE 

Code de la propriete intellectuelle - Livre VI 



N° 11235*02 



DEP/(RTEMENT DES BREVETS 

26 bis, rue de Saint Petersbourg 
75800 Paris Cedex 08 

Telephone : 01 53 04 53 04 Telecopie : 01 42 93 59 30 



DESIGNATION D'INVENTEUR(S) Page N** 

(Si le demandeur n'est pas I'inventeur ou Tunique inventeur) 



Cet imprim6 est a remplir lisiblement a Tencre noire 



DB 113W/260899 



Vos references pour ce dossier 

(facuitatif) 


JC/CLN 


N° D'ENREGISTREMENT NATIONAL 





TITRE DE L'iNVENTION (200 caractdres ou espaces maximum) 

METHODE DE MAILLAGE AUTOMATIQUE DE CONDUITES PERMETTANT L'IMPLEMENTATION DE CODES DE 
MODELISATION DE FLUX 



LE(S) DEMANDEUR(S) : 

INSTITUT FRANCAIS DU PETROLE 



DESIGNE(NT) EN TANT QU'INVENTEUR(S) : (Indiquez en haut a droite «Page N° 1/1» S'il y a plus de trois inventeurs, 
utilisez un formulaire identique et numerotez chaque page en indiquant le nombre total de pages). 



Norn 


DURET 


Prenoms 


Emmanuel 


Adresse 


Rue 


8 rue Isabey 




Code postal et vilte 


92500 Rueil Malmaison 


Societe d'appartenance (facultatij) 




Nom 


FAILLE 


Prenoms 


Isabella 


Adresse 


Rue 


1 Le Val Eglantine 




Code postal et ville 


78240 Carriere su Seine 


Societe d'aooartenance (facuitatif) 




Nom 


HEINTZE 


Prenoms 


Eric 


Adresse 


Rue 


3 rue du Progres 




Code postal et ville 


92190 Meudon 


Societe d'appartenance (facuitatif) 




DATE ET SIGNATURE(S) 
DU (DES) DEMANDEUR(S) 
OU DU MANDATAIRE 
(Nom et qualite du signataire) 

Alfred ELMALEH, >0«=^ 
Chef du Departement Brevets^^*^''^!^ . 
N° 422-5/PP.253 \^ 





La loi n°78-17 du 6 janvier 1978/^ative a I'informatique, aux fichiers et aux libertes s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de I'lNPI. 





DOCUMENT COMPORTANT DES MODIFICATIONS 



PAGE(S) OE LA DESCRIPTION OU DES REVENDICATIONS 
OU PLANCHE(S) DE DESSIN 


R.M/ 


DATE 
OELA 
CORRESPONDANCE 


TAMPON DATEUR 
OU 

CORRECTEUR 


Modifiee(4) 


Supprimee(s) 


Ajoutefl(s) 










ohhoioo 


































































































1 



Un changement apporte a la redaction des revendications d'origine. sauf si ceiunci d6coule des dispositions de rarticle R.612.36 
du code de la Propri6t6 Inteilectuelle, est signal^ par la mention cR.M.* (revendicaUons modifees). 




# 



1 



La presente invention conceme une methode de maillage automatique de 
conduites permettant T implementation de codes de modelisation de fluides 
achemines par ces conduites. 

5 La methode selon 1' invention trouve des applications dans de nombreux 

domaines. On peut notamment I'utiliser dans le domaine de la production 
d'hydrocarbures pour Timplementation de codes de simulation d'ecoulements 
polyphasiques dans des conduites petrolieres depuis des sites de production jusqu'a 
des sites de destination. 

10 Le maillage obtenu par la methode peut etre utilise notamment pour 

r implementation du code de modelisation TACITE (marque deposee) destine a 
simuler des ecoulements permanents ou transitoires d'hydrocarbures dans des 
conduites. Differents algorithmes permettant de conduire la simulation des 
ecoulements suivant le code TACITE ont fait I'objet des brevets US 5 550 761, 

15 FR 2 756 044 et FR 2 756 045 (US 5 960 187). 

Les modes d'ecoulement de fluides polyphasiques dans des tubes sont 
extremement varies et complexes. Les ecoulements diphasiques, par exemple, 
peuvent etre stratifies, la phase liquide s'ecoulant dans la partie inferieure de la 
conduite, ou intermittents avec une succession de bouchons liquides et gazeux, ou 

20 bien encore disperses, le liquide etant entraine sous forme de fines gouttelettes. Le 
mode d'ecoulement varie notamment avec I'inclinaison des conduites par rapport a 
I'horizontale et il depend du debit de la phase gazeuse, de la temperature etc. Le 
glissement entre les phases qui varie selon que Ton considere les trongons de 
conduite ascendants ou descendants, entraine des variations de pression sans qu'il y 

25 ait toujours compensation. Les caracteristiques du reseau d'ecoulement 
(dimensionnement, pression, debit de gaz, etc.) doivent etre determinees avec soin. 



Le code de simulation TACITE prend en compte un certain nombre de 
parametres influant directement sur la physique du probleme a trailer. Parmi ceux-ci, 
on trouve les proprietes des fluides et des modes d'ecoulement, les variations 
topographiques (variations de longueurs, d'inclinaison, de diametre etc.), la rugosite 
eventuelle des conduites, ses proprietes thermiques (nombre de couches d'isolant et 
leur nature) ou encore la disposition d'equipements le long de la conduite (pompes, 
injecteurs, separateurs, etc.) qui entrainent des changements de flux physique. 

Etat de la technique 

Le maillage d'un domaine physique est une etape primordiale lors d'une 
simulation numerique, De sa qualite dependent la validite des resultats et les temps 
de calcul. II est done fondamental de foumir au code un maillage correct avant de 
commencer une simulation. On juge generalement la qualitd d'un maillage a sa 
capacite a bien decrire les phenomenes physiques sans que la simulation ne prenne 
trop de temps, si bien qu'il existe toujours un maillage optimal pour chaque probleme 
etudie. Un maillage inadapte peut conduire, dans I'execution du schema numerique 
qui regit la simulation, a des erreurs difficilement detectables, du moins au premier 
abord, voire meme rendre le calcul impossible et stopper Texecution du code, s'il est 
par trop aberrant. Les utilisateurs des codes n'ont pas forcement une experience 
suffisante en analyse numerique pour la realisation d'un maillage correct, susceptible 
de rendre reellement compte des phenomenes physiques que Ton veut etudier. 

La topographic d'une conduite cylindrique peut etre assimilee a une 
succession de segments de droites reliant des points successifs. En coordonnees 
cartesiennes, deux points successifs de la conduite sur ses portions verticales 
(montantes ou descendantes) peuvent avoir meme abscisse (courbe A sur la Fig.l). II 
est bien preferable de ce fait de representer Televation de chaque point en fonction de 
son abscisse curviligne le long de la conduite. Avec ce mode de representation, des 
points successifs de la conduite d'elevations differentes ont forcement deux abscisses 
curvilignes distinctes et la pente des tron^ons de conduite est au plus egale a 45° par 
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rapport a I'horizontale (cas des trongons montants ou descendants absolument 
verticaux courbe B de la Fig.l). A une abscisse correspond toujours une ordonnee et 
une seule. 

Avec un peu de sens physique, certaines erreurs de maillage peuvent etre 
5 evitees. Aux endroits de la conduite susceptibles de connaitre de grandes variations 
des parametres physiques si on est en mesure de les prevoir, on peut imposer un 
maillage plus fin qu'ailleurs. Ainsi, on effectue moins de calculs a chaque pas de 
temps, tout en conservant la finesse souhaitee aux endroits importants. Mais, dans 
I'optique d'obtenir une solution continue, le passage d'une maille fine a une maille 
10 plus grossiere doit, lui aussi, etre continu. 

La Fig. 2a montre par exemple une portion de conduite en W d'une longueur 
de 2 km comprenant quatre tron§ons de 500 m. Si on discretise une telle conduite par 
des mailles a pas de 40 m constant du debut a la fin, on laisse de cote les points 
importants du trace a 500 m et 1500 m. La simulation ne permettra pas de rendre 
15 compte correctement de 1' accumulation de liquide a ces point bas de la topographic. 
Plus important encore, le calcul est fausse par le fait qu'on ait remplace les angles du 
W par des segments de droites horizontaux (Fig. 2b). Les phenomenes physiques ne 
seront pas ceux recherches. 

La methode selon Tinvention permet d'obtenir automatiquement un maillage 
20 ou discretisation d'une conduite tenant compte au mieux de la topographie et des 
parametres physiques affectant la physique des ecoulements, soumis aux contraintes 
suiv antes ; 

■ 1 - Assurer la convergence du calcul ; 

■ 2 - Representer au mieux les accumulations importantes de liquide dans les 
25 points bas de la conduite ; 

■ 3 - Placer les equipements sur un bord de maille ; 

■ 4 - Imposer a deux mailles consecutives d'avoir le meme ordre de longueur ; 

■ 5 - Respecter la longueur totale de la conduite ; 
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■ 6 - Limiter le nombre de mailles au minimum possible en respectant les 
contraintes precedentes, pour ne pas trop penaliser la simulation avec le temps de 
calcul. 

Le respect des six contraintes precedentes n'est pas une tache aisee, mais il est 
5 le passage oblige pour eviter de mailler la conduite etudiee de fagon homogene, sans 
se soucier de la physique du probleme, comme le font la plupart des mailleurs 
automatiques. 

Pour limiter le nombre de mailles, il faut chercher, si cela est possible, a 
simplifier la topographie pour ne retenir que les seules zones de la conduite ou se 
10 trouvent les variations significatives du profil susceptibles d'intervenir de fa?on 
significative sur les phenomenes physiques. 

La methode selon i'invention 

La methode selon I'invention permet de mailler automatiquement une 
conduite ID presentant sur toute sa longueur une topographie ou profil quelconque 
15 dans le but de faciliter 1' implementation de codes de modelisation de flux. Le 
maillage obtenu par la methode presente une repartition de mailles de dimensions 
variables, appropriee pour tenir compte au mieux de la physique des ecoulements. 

La methode est caracterisee en ce que, ayant defini une taille minimale de 
maille et une taille maximale de maille, on subdivise la conduite en trongons 
20 delimites par des coudes, on positionne une maille de taille minimale de part et 
d' autre de chaque coude, on positionne des grandes mailles de taille au plus egale a 
la taille maximale dans la portion centrale de chaque tron9on, et on repartit des 
mailles de tailles croissantes ou decroissantes sur les portions intermediaires de 
chaque tron^on entre chaque maille de taille minimale et la portion centrale. 

25 La repartition des mailles de tailles croissantes ou decroissantes sur les 

portions de chaque tron9on intermediaire entre chaque maille de taille minimale et la 
portion centrale, est obtenue par exemple en determinant les points d' intersection 
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avec chaque tron^on de conduite, d'un faisceau de droites concourant en un point et 
formant entre elles un angle constant. 

On detemnine par exemple la position du sommet du faisceau de droites sur 
un axe passant par un coude de la conduite et perpendiculaire a chaque trongon, a 
5 une distance de celui-ci qui est fonction de la taille des mailles extremes de chaque 
portion intermediaire et de leur ecart. 

Le positionnement automatique des mailles avec des mailles plus petites au 
voisinage des extremites de chaque trongon permet d'apporter un plus grand soin a la 
modelisation des phenomenes dans les portions de conduite presentant des 
10 changements de direction (inflechissement ou coude). 

La methode selon 1' invention comporte de preference une simplification 
prealable de la topographic de la conduite de fa^on que le nombre total de mailles du 
maillage de la conduite permettent une modelisation realiste de la physique des 
phenomenes dans un temps fixe. 

15 Suivant un premier mode d' implementation, la methode comporte une 

representation de la conduite sous la forme d'un graphe reliant I'abscisse curviligne 
et la variation de niveau, et une simplification du nombre de tron^ons a) en affectant 
a chaque point entre deux tron9ons successifs un poids prenant en compte la 
longueur des tron§ons et leurs pentes respectives, b) en selectionnant parmi les points 

20 ranges par ordre de poids croissant ou decroissant, ceux dont le poids est le plus 
eleve, la topographic simplifiee etant celle du graphe passant par les points 
selectionnes. 

La selection des points de la conduite dont le poids est le plus eleve est 
obtenue par exemple en reperant dans le rangement de points une discontinuite de 
25 poids superieure a un certain seuil fixe. 

Suivant un autre mode d'implementation, la methode comporte une 
representation de la conduite sous la forme d'un graphe reliant I'abscisse curviligne 



et la variation de niveau, et une simplification du nombre de tron9ons par a) la 
formation du spectre de frequence de la courbe representative de la topographie de la 
conduite, b) Tattenuation des plus hautes frequences du spectre traduisant les plus 
petites variations de topographie et c) la reconstruction d'une topographie simplifiee 
correspondant au spectre de frequence rectifie. 

La selection est effectuee par exemple par a) un echantillonnage de la courbe 
representative de la topographie de la conduite avec un pas d' echantillonnage choisi 
pour que le plus petit tron^on de la conduite contienne au moins deux pas 
d'echantillonnage, b) une determination du spectre de frequence de la courbe 
echantillonnee par application, c) une correction du spectre par filtrage passe-bas 
dont la frequence de coupure est choisie en fonction d'un nombre de mailles 
maximal fixe pour subdiviser la conduite, et, d) une determination de la topographie 
correspondant au spectre de frequence rectifie. 

Les deux modes de simplification automatique precedents peuvent etre 
appliques independamment I'un de Tautre ou bien successivement, le second mode 
etant applique de preference quand le premier mode ne permet pas de simplification 
notable de la topographie. 

Presentation sommaire des figures 

D'autres caracteristiques et avantages de la methode selon I'invention, 
apparaitront a la lecture de la description ci-apres d'exemples non limitatifs de 
realisation, en se referant aux dessins annexes ou : 

-la Fig.l montre deux modes de representation schematique de la topographie 
d'une conduite ; 

les Fig. 2a, 2b montrent respectivement la topographie schematique d'une 
conduite en forme de W en coordonnees curvilignes, et une partie agrandie de 
cette meme topographie ; 



la Fig. 3 montre un mode d' attribution de poids a un point de la topographic 
d'une conduite ; 

la Fig. 4 montre un example de spectre de poids ; 

la Fig. 5 montre un exemple d' arrangement des points par paliers de poids 
decroissants, permettant de reperer la position d'un seuil et de simplifier la 
topographic de la conduite ; 

la Fig. 6 montre un exemple de topographic d'une conduite marine comportant 
un « riser » a ses extremites ; 

la Fig. 7 montre la topographic simplifiee de la meme conduite, que Ton obtient 
par selection des poids ; 

la Fig. 8 montre que, sans les « risers » terminaux, la forme generale de la meme 
conduite est plus difficile a degager ; 

la Fig. 9 montre un spectre de frequence typique d'une conduite ; 

la Fig. 10 montre un exemple de tron9on de conduite avec une distribution de 
mailles de differentes tailles, les plus petites Ml etant positionnees aux coudes, 
les plus grandes, M2 etant placees dans le tiers central, les mailles intermediaires 
M3 etant interposees et resultant d'une interpolation I entre les unes et les autres ; 

la Fig. 11 montre un mode de formation de mailles de tailles croissantes ; 

la Fig. 12 illustre le mode de decoupage angulaire d'une portion intermediaire sur 
un trongon de conduite ; et 

la Fig. 13 montre le maillage obtenu par la mise en ocuvre de la methode, sur une 
conduite sous-marine de 90km de longueur. 
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DESCRIPTION DETAILLEE 

I) Simplification de la topographie d'une conduite 

On n'a generalement aucune difficulte a faire ressortir au premier coup d'ceil 
la forme globale d'un profil quelconque. La methode selon rinvention permet, par 
5 des criteres purement mathematiques, le reperage automatique de la configuration 
d'une conduite, en se basant sur une analyse spectrale de la courbe representative des 
variations du profil. Parmi tous les spectres que Ton peut associer a une topographie 
donnee, on en recherche un qui permette de distinguer les portions du profil a 
simplifier et les portions du profil importantes. 

10 I-l) Premier mode de simplification 

Dans une topographie, les seuls criteres faisant en sorte qu'un point est 
simplifiable par rapport a un autre, ne peuvent etre que les longueurs des sections qui 
I'entourent et la difference angulaire les separant (Fig.3). Lorsqu'on construit les deux 
« spectres » (Indices du trongon) - (Longueurs des trongons) et (Abscisse curviligne 
15 des points) - (Difference angulaire des tron9ons entrant et sortant), on se rend compte 
qu'ils presentent des differences notables d'ordres de grandeur mais egalement que 
ces deux spectres sont independants, si bien qu'en simplifiant des points negligeables 
dans I'un, il se peut qu'on ait efface des points importants dans I'autre. 

Pour regrouper ces deux spectres en un seul, on va affecter a chaque point 
20 topographique un poids prenant en compte les longueurs de trongons, et les 
differences angulaires qui les separent. On utilise par exemple la ponderation 



suivante : 



Poids — 



(Pt -/*,)■ ou Li et Lj sont les longueurs des tron§ons, et 




et P^=— sont les pentes. Ainsi, a longueurs egales, on simplifiera les 
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tron^ons separes par la plus petite difference de pente. Et a angles egaux, on 
simplifiera les longueurs les plus courtes. 

Construction du spectre 

Le spectre (Abscisse curviligne - Poids) presente dans la majorite des cas, une 
5 succession de pics de toutes tailles. Ces spectres tel celui de la Fig.4 ne sont 
generalement pas analysables directement. Dans ces conditions, la technique que Ton 
utilise ici consiste a classer les poids par ordre croissant ou decroissant et a leur 
affecter a chacun Tindice correspondant de classement par poids de 1 a N. De 
preference, on utilise une representation (Log Poids- Indice) qui fait mieux ressortir 
10 les ordres de grandeur car un saut de n sur un tel spectre signifie un rapport de 10" 
sur les poids. Tous les poids ayant le meme ordre de grandeur sont classes sur des 
paliers plus ou moins horizontaux. Deux poids d'ordres de grandeur differents seront 
separes par un segment de droite vertical. On aboutit a un . spectre en cascade, 
permettant de lire aisement les differents ordres de grandeurs presents dans la 
15 topographie. Sur I'exemple de la Fig.5 par exemple, le spectre logarithmique contient 
deux paliers bien distincts separes par un segment vertical. 

On cherche le premier triplet de points consecutifs du spectre, defini par 
exemple par un seuil AP fixe sur I'echelle logarithmique ( AP =1 par exemple) entre 
le deuxieme et le troisieme, qui suit un saut inferieur a AP entre le premier et le 

20 deuxieme. Les deux premiers points sont du meme ordre de grandeur. Tous les points 
suivants sont d'un ordre de grandeur negligeable par rapport aux deux premiers . De 
cette fa9on, on s'assure que tous les poids a droite du triplet en question, seroni au 
moins 10 fois inferieurs au poids du deuxieme et done negligeables par rapport aux 
points en amont. Sur la table de correspondance (indice de poids-abscisse 

25 curvihgne), on selectionne les points d'abscisse curviligne correspondant aux poids 
les plus forts selectionnes. La topographie simplifiee sera la ligne passant par ces 
points. 
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Sur I'exemple de topographic de la Fig. 6, on y distingue trois parties 
distinctes. Elle commence par un « riser » de 3 kms, suivi d'une partie horizontale en 
dents de scie sur 20 kms qui se termine par un autre « riser » de 200 m lui aussi en 
dents de scie. Son spectre est celui de la Fig.5. Le premier triplet qui obeit au critere 
5 de seuillage est forme par les points 4, 5 et 6. Le seuil de simplification est le point 
d'indice 6. Un saut superieur a 2 dans Techelle logarithmique separe les paliers 
horizontaux de part et d'autre des points 5 et 6. Ainsi, on s*assure que les points a 
gauche de Tindice 5 ont des poids au moins 100 fois superieurs a ceux qui se trouvent 
a droite de I'indice 6. 

10 Dans cet exemple, on simplifie la topographic en ne gardant que les points 

d'abscisse curviligne correspondant aux poids superieurs ou egaux a celui du point 6. 
On obtient la topographic simplifiee de la Fig. 7. La forme globale est conservee. 
Toutes les petites variations en dents de scie sur la partie horizontale de 20 kms ont 
ete supprimees. Le nombres de points est passe de 43 initialement (Fig. 6) a 6, soit 

15 une reduction d'un facteur 7. Ce cas sc prete particulierement bien au seuillage, 
puisque les differcnts ordres de grandeurs sont visiblcs sur la topographic initialc. 

Le premier mode de simplification qui vient d'etre decrit est facile a mettre en 
oeuvre et fonde sur des algorithmes relativement simples executables rapidement. II 
est adapte aux topographies possedant plusieurs ordres de grandeurs, comme la 
20 topographic precedente qui a pu etre considerablement simplifiee car elle contenait 
des points ayant des poids negligeablcs par rapport a d'autres. 

Le probleme est tout autre si on ne s'interesse qu'a la partie centrale de cette 
topographic, en retirant les risers aux extremites car dans ce cas, comme le montre la 
Fig.8, la forme generalc de la conduite est plus difficile a degager. La simplification 
25 de cette topographic par une droite reliant le point de depart et celui d'arrivee, n'est 
plus possible. Son spectre est exactement le meme que celui de la topographic 
initiate, a ceci pres qu'il commence au point 6. Aucun seuil n'est present dans cette 
partie du spectre, les points ont tous le meme ordre de grandeur. Et meme si le poids 
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le plus eleve vaut plus de 100 fois le plus faible, on passe de Tun a I'autre de fagon 
continue. 

1-2) Deuxieme mode de simpliflcation 

Pour les topographies avec des points presentant le meme ordre de grandeur, 
5 qui ne peuvent pas etre traitees par la methode de seuillage precedente, on procede a 
un filtrage spectral. Les petites variations dans le profil de la conduite se traduisent 
par des hautes frequences dans le spectre de Fourier de la fonction representative de 
la topographic. En coupant ou en attenuant les frequences les plus hautes de son 
spectre de frequence, il est possible de simplifier la topographic. 

10 A cet effet, on echantillonne la fonction topographique et Ton determine son 

spectre par la methode dite FFT (Fast Fourier Transform). Le pas d*echantillonnage 
doit etre assez petit pour bien rendre compte de toutes les plages de frequence en 
evitant le repliement du spectre. Pour ce faire, le nombre de points d'echantillonnage 
est choisi de telle fa^on que le plus petit trongon de conduite contienne au moins 

15 deux subdivisions pour s'assurer que la transformee de Fourier agira sur toutes les 
parties de la conduite, meme les plus insignifiantes.L'attenuation des hautes 
frequences doit bien entendu etre effectuee avec discemement et ajustee de fagon que 
la fonction topographique reconstituee reste representative de la fonction initiale. 

La methode de filtrage la plus simple consiste par exemple a appliquer un 
20 seuil, tous les coefficients de Fourier inferieurs a ce seuil etant elimines (coefficients 
inferieurs a 40 par exemple sur Texemple de la Fig.9. Seule I'information contenue 
dans les frequences inferieures a ce seuil est conservee. On reconstitue par 
transformee inverse la typographic simplifiee correspondante. 

En se fixant une frequence de coupure, on fixe ainsi le nombre maximal 
25 d'oscillations du signal reconstitue. Si Ton ne garde que les dix premieres frequences, 
la fonction reconstituee va suivre la forme generale de la conduite, avec un 
maximum de vingt extrema. 
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II) Selection des tallies de maille sur chaque tron^on de conduite 

Principe 

Le principe du maillage va consister a mailler independamment les trongons 
de conduite entre deux bords imposes. Comme I'interet d'un maillage bien fait est de 

5 pouvoir observer correctement les accumulations de liquide dans les coudes, il est 
preferable que le maillage soit affine aux points de la topographie susceptibles de 
voir s'accumuler du liquide ou du gaz, C'est pourquoi, on va s'arranger pour placer 
une maille courte avant et apres chaque coude et des plus grandes entre les coudes. 
Par contre, il n'est pas utile de mailler finement les parties intermediaires des 

10 tron§ons entre les coudes. 

La topographie de la conduite ayant ete au prealable (quand c'etait 
necessaire) simplifiee et reduite a un certain nombre de tron^ons, on fixe pour les 
mailles une taille minimale et une taille maximale. On isole alors tout d'abord les 
bords de chacun d'eux (leur entree, leur sortie) par des cellules de petite taille et 
15 ensuite on insere des bords de mailles sur leur partie centrale qui est plus longue. II 
n'est generalement pas utile d'af finer le maillage en entree et en sortie en dehors des 
portions aux extremites de chaque tron^on et Ton pent done inserer sur une bonne 
part de la longueur de chaque trongon (par exemple sur les 2/3 de la longueur) de la 
taille maximale que Ton s'est fixe. 

20 Comme distribution, on peut choisir par exemple que la taille des mailles 

apres celle qui suit un coude, augmente graduellement sur un tiers de la longueur du 
trongon, reste constante sur le tiers suivant pour enfin diminuer graduellement sur le 
dernier tiers avant la maille courte finale comme le montre la Fig. 10. 

Definition des longueurs de cellules minimale et maximale 

25 On va definir deux longueurs de cellules, une longueur minimale, qui va 

servir a isoler les bords de mailles imposes par des cellules de petites tailles, et une 
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longueur maximale, qui va servir a mailler les milieux des trongons compris entre 
deux cellules couites. 

Toutes les mailles qu*on insere apres ces deux etapes sont deduites des 
cellules initiales par interpolation entre une cellule courte et une cellule longue. Elles 
5 ont done des tailles intermediaires. Cette propriete est interessante. Elle indique que 
le nombre de mailles total sera forcement compris entre le nombre qu'on aurait 
obtenu en maillant de fa9on homogene avec la longueur minimale, et le nombre 
obtenu de la meme fagon mais avec la longueur maximale. Ainsi, on peut controler le 
nombre de mailles total a partir des tailles minimale et maximale. 

10 Une des contraintes du maillage automatique reside dans le nombre de 

mailles total. Celui-ci doit engendrer un temps de simulation le plus court possible, 
tout en permettant une bonne visualisation des phenomenes physiques^ L'experience 
montre d'une part qu'une discretisation de moins de 40 mailles ne permet pas une 
bonne description physique des problemes. D'autre part, les maillages de plus de 150 

15 mailles engendrent des simulations trop longues. Le maillage par defaut doit done 
etre assez souple et comporter entre 40 et 100 mailles. 

Un nombre de mailles aussi faible ne convient pas pour tous les cas. Le 
nombre de mailles ideal pour un cas precis depend de plusieurs facteurs, entrant en 
compte dans le schema numerique. A topographie egale par exemple, un cas 
20 comportant un grand nombre de changements de section necessitera un maillage plus 
fin. La methode selon 1' invention laisse une grande latitude a Tutilisateur de choisir 
le nombre total de mailles qui lui convient. 

A partir de ce nombre N, le code calcule les longueurs minimales Min et 
maximale Max de la fagon suivante : 

L 

25 Min = 
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Max = ^ 



N-P 



Le parametre P permet de reduire Tecart entre les longueurs minimale et 
maximale de fa^on a rendre progressivement le maillage homogene pour ]es grands 
nombres de mailles. 

5 On definit ce parametre par exemple comme suit. Pour un nombre de mailles 

demande inferieur ou egal a 60 par exemple, on le fixe par exemple a 60. C'est le 
maillage par defaut. Le parametre vaut 40. La plus petite maille vaudra L/lOO et la 
plus grande L/20. Le nombre de mailles total sera compris entre 20 et 100. 

Un nombre de mailles superieur ou egal a 150 signifie que la modelisation a 
10 traiter est certainement plus delicate. Le but est alors de construire un maillage 
homogene. Pour cela, les tallies minimales et maximales doivent etre'proches Tune 
de I'autre. On va done fixer le parametre a 10. Le nombre de mailles total sera alors 

compris entre — — — et — — — .On obtient le nombre de mailles demande, a 20 
A^ + IO A^-10 

mailles pres, sur plus de 150. 

15 Pour faire en sorte que le maillage s'homogeneise progressivement entre 60 et 

150 mailles, on va calculer le parametre par interpolation lineaire entre les deux 
domaines, ce qui s'exprime de la fa^on suivante : 

P=40 si N<60 

P=-1n + 60 si 60<N<150et 
3 

20 P=10 si N>150 



Ce parametre etant determine, il est possible d'isoler les bords imposes par 
des cellules courtes et de discretiser les milieux des trongons par des cellules 
longues. 
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II ne reste plus qu*a trouver un moyen pour passer graduellement d'une cellule 
courte a une longue. On connait les longueurs des trois cellules, et on cherche a 
inserer des bords de mailles sur la partie centrale. Les tallies des mailles ainsi creees 
doivent etre comprises entre celles des cellules extremes. Partant de la plus petite, la 
maille suivante doit toujours etre plus longue que la precedente, mais plus courte que 
la suivante. 

Dans le cas general, il n'existe aucun couple (/,n) e (/?, A^) tel que : 

- la taille d^une maille se deduise de celle de la precedente en la multipliant par un 
facteur f . 

la somme des n longueurs ainsi creees soit egale a (L1+L2) 

la taille de la demiere cellule puisse s'exprimer sous la forme /"""VLi f . 

II en va de meme pour une eventuelle interpolation lirieaire entre les deux 
mailles. Le fait de connaitre les trois longueurs impose une surabondance de donnees 
par rapport aux inconnues. II est alors impossible de satisfaire toutes les contraintes. 

Pour pallier cette difficulte, on propose une methode de type geometrique ou 
Ton utilise la propriete que les segments LI, L2, L3, L4 decoupes sur un axe par les 
droites d'un faisceau regulier (d'ecartement angulaire a constant les unes par rapport 
aux autres) dont le sonrmiet est en dehors de cet axe, varient progressivement 
(Fig. 11). 

On considere (Pig. 12) un tron9on de conduite commengant par une petite 
maille (0, xl) de longueur LI et se terminant par une maille (x2, x3) de longueur 
L3 > LI. On peut montrer qu'il existe un point sur une perpendiculaire au trongon de 
conduite a I'abscisse 0 tel que les mailles de longueurs LI et L3 soient vues de ce 
point sous le meme angle a. L'ordonnee y de ce sommet est donne par la relation : 
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ou L2 est la longueur du segment (x 1 , x2). 

II s'agit ensuite de couper Tangle P en N parties egales, N etant egal a la 

division entiere de B par a, soit N= £"(—). Chacun des N angles decoupant P est 

a 

toujours superieur ou egal a a. 

Le principe utilise pour realiser T insertion des bords de mailles est a la fois 
simple et souple. D permet, par le biais d*un unique parametre, de creer un maillage 
soit uniforme, soit heterogene affine aux endroits importants. 
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REVENDICATIONS 

1) Methode de maillage automatique de conduites permettant 
r implementation de codes de modelisation de fluides achemines par ces conduites, 

5 caracterisee en ce que, ayant defini une taille minimale de maille et une taille 
maximale de maille, on subdivise la conduite en tron^ons delimites par des coudes, 
on positionne une maille de taille minimale de part et d' autre de chaque coude, on 
positionne des grandes mailles de taille au plus egale a la taille maximale dans la 
portion centrale de chaque tron^on, et on repartit des mailles de tallies croissantes ou 

10 decroissantes sur les portions intermediaires de chaque trongon entre chaque maille 
de taille minimale et la portion centrale. 

2) Mdthode selon la revendication 1, caracterisee en ce que Ton repartit des 
mailles de tallies croissantes ou decroissantes sur les portions de chaque trongon 
intermediaire entre chaque maille de taille minimale et la portion centrale en 

15 determinant les points d' intersection avec chaque tron9on de conduite, d'un faisceau 
de droites concourant en un point et formant entre elles un angle constant. 

3) Methode selon la revendication 1, caracterisee en ce que Ton determine la 
position du sommet du faisceau de droites sur un axe passant par un coude de la 
conduite et perpendiculaire a chaque tronqon, a une distance (y) de celui-ci qui est 

20 fonction de la taille (L'l, L'3) des mailles extremes de chaque portion intermediaire 
et de leur ecart (L'2). 

4) Methode de maillage selon Tune des revendications 1 a 3, caractensde en 
ce qu'elle comporte une simplification prealable de la typographique de la conduite. 

5) Methode de maillage selon la revendication 4, caracterisee en ce qu'elle 
25 comporte une representation de la conduite sous la forme d'un graphe reliant 

I'abscisse curviligne et la variation de niveau, et une simplification du nombre de 
tron9ons en affectant a chaque point entre deux trongons successifs un poids prenant 
en compte la longueur (LI, L2) des tron9ons et leurs pentes respectives (PI, P2) et en 
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selectionnant parmi les points ranges par ordre de poids croissant ou decroissant, 
ceux dont le poids est le plus eleve. 

6) Methode de maillage selon la revendication 5, caracterisee en ce Ton 
selectionne les points de la conduite dont le poids est le plus elev6 en reperant dans le 

5 rangement de points une discontinuite de poids superieure a un certain seuil fixd 
(AP). 

7) Methode de maillage selon la revendication 5, caracterisee en ce qu'elle 
comporte une representation de la conduite sous la forme d'un graphe reliant 
I'abscisse curviligne et la variation de niveau, et une simplification du nombre de 

10 trongons par la formation du spectre de frequence de la courbe representative de la 
topographie de la conduite, T attenuation des plus hautes frequences du spectre 
traduisant les plus petites variations de topographie et la reconstruction d'une 
topographie simplifiee correspondant au spectre de frequence rectifie. 

8) Methode de maillage selon la revendication 7, caracterisee en ce qu'elle 
15 comporte un echantillonnage de la courbe representative de la topographie de la 

conduite avec un pas d'echantillonnage choisi pour que le plus petit trongon de la 
conduite contienne au moins deux pas d'echantillonnage, une determination du 
spectre de frequence de la courbe echantillonnee par application, une correction du 
spectre par filtrage passe-bas dont la frequence de coupure est choisie en fonction 
20 d'un nombre de mailles maximal fixe pour subdiviser la conduite, et la determination 
de la topographie correspondant au spectre de frequence rectifie. 
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REVENDICATIONS 

1) Methode de maillage automatique de conduites permettant 
r implementation de codes de modelisation de fluides acheminds par ces conduites, 

5 caracteris^e en ce que, ayant defmi une taille minimale de maille et une taille 
maximale de maille, on subdivise la conduite en trongons delimites par des coudes, 
on positionne une maille de taille minimale de part et d' autre de chaque coude, on 
positionne des grandes mailles de taille au plus egale a la taille maximale dans la 
portion centrale de chaque trongon, et on repartit des mailles de tailles croissantes ou 

10 decroissantes sur les portions intermediaires de chaque tron^on entre chaque maille 
de taille minimale et la portion centrale. 

2) Methode selon la revendication 1, caracterisee en ce que Ton repartit des 
mailles de tailles croissantes ou decroissantes sur les portions de chaque tron§on 
intermediaire entre chaque maille de taille minimale et la portion centrale en 

15 determinant les points d' intersection avec chaque tron^on de conduite, d'un faisceau 
de droites concourant en un point et formant entre elles un angle constant. 

3) Methode selon la revendication 1, caracterisee en ce que Ton determine la 
position du sommet du faisceau de droites sur un axe passant par un coude de la 
conduite et perpendiculaire a chaque tron^on, a une distance (y) de celui-ci qui est 

20 fonction de la taille (LI, L3) des mailles extremes de chaque portion intermediaire et 
de leur ecart (L2). 

4) Methode de maillage selon Tune des revendications 1 a 3, caracterisee en 
ce qu'elle comporte une simplification prealable de la topographic de la conduite. 

5) Methode de maillage selon la revendication 4, caracterisee en ce qu'elle 
25 comporte une representation de la conduite sous la forme d'un graphe reliant 

I'abscisse curviligne et la variation de niveau, et une simplification du nombre de 
tron^ons en affectant a chaque point entre deux tron§ons successifs un poids prenant 
en compte la longueur (LI, L2) des trongons et leurs pentes respectives (PI, P2) et en 
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